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Recent progress on isospin breaking corrections and 
their impact on the muon g — 2 value* 
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Abstract We describe some recent results on isospin breaking corrections which are of relevance for predictions 
of the leading order hadronic contribution to the muon anomalous magnetic moment a^ ad ' LO when using r 
lepton data. When these corrections are applied to the new combined data on the tt^tt spectral function, the 
prediction for aJj ad,LO based on r lepton data gets closer to the one obtained using e + e~ data. 
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1 Introduction 

Currently, the accuracy of the Standard Model 
prediction of the muon anomalous magnetic moment 
a M = (g — 2) /2 is limited by the uncertainties of 
hadronic contributions [l|, |2| . The dominant term in 
the leading order hadronic contribution a h ad > LO anc j 
an important part of its associated uncertainty is pro- 
vided by the 7T7t spectral function, which can be mea- 
sured in e + e~ annihilations and in t lepton decays 
(more details about their current status are given 
in the accompanying contribution by Michel Davier 
[2]). Owing to the isotopic properties of the elec- 
tromagnetic and AS = weak vector currents, the 
so-called Conserved Vector Current (CVC) hypothe- 
sis, the spectral functions themselves and their con- 
tributions to flh ad ' LO must be the same after isospin 
breaking (IB) corrections are appropriately applied 
to input data B- 

In recent years, a comparison of e + e~ and r based 
measurements of the 7T7t spectral functions in the 
timclike region, have shown large discrepancies for 
center of mass energies above the p(770) resonance 
peak, beyond the size expected for IB corrections 
0, B- The predictions for a h ad - LO based on these 
two sets of data have been in disagreement by more 
than 2a 

BB- 

Moreover, the branching fraction for 
t — > ttttv predicted from e + e~ data corrected by IB 
effects was underestimated by more than 4tr with re- 
spect to the average of direct measurements B- 



Given this 'e + e~ vs r discrepancy'*, it is believed 
that t decay data does not provide at present a re- 
liable determination of a^^ (currently, other use- 
ful contributions from r decay data involve only the 
27T and 4-7T channels 0]). Even though unidentified 
errors may be affecting e + e~ and/or t decay data, 
understanding IB corrections becomes crucial to gain 
confidence about this important 27r contribution and, 
when consistency is achieved, to have a more precise 
prediction of a^ad.LO f rom combined results. 

In this contribution we summarize some recent re- 
sults on the isospin breaking corrections that are rel- 
evant for understanding such discrepancies. As it is 
discussed in @, 0], their application to the evalua- 
tion of a^ ad,LO from r data leads to a value B that 
is closer to e + e"-based calculations. The prediction 
of the r — > ttttv branching fraction based on the IB 
corrected e + e~ data also shows a reduced discrepancy 
with respect to results of direct measurements [2J, |7J . 
A preliminary version of this work appeared in rcf. 
8] . In this new version we include a discussion of the 
complete set of IB corrections and we address some 
points that were unclear in some of our previous re- 
ports. 

2 IB corrections to the a^ adLO disper- 
sion integral from r data 

The leading order hadronic contribution a h ad ' LO 
can be evaluated by using a combination of experi- 
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'Another reason is that e+e data is more directly related to aJj ad ' LO through the dispersion integral than r decay data. 
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mental data and perturbative QCD for the hadronic 
vacuum polarization (HVP) function of the photon. 
At low energies, where QCD does not provide a reli- 
able calculation of Green functions, the HVP can be 
constructed as a sum over exclusive hadronic chan- 
nels measured in e + e" annihilation. The dispersion 
integral relating each exclusive e + e~ — > X° channel 



to a^' LO is: 



(1) 



4 /n 'i 



where R^l (s) is the ratio of hadronic X° to pointlike 
[i + \i~ bare cross sections [H, E| in e + e~ annihilation 
at a center of mass energy yfs. The behavior of the 
QED kernel K(s) ~ l/s [9j, enhances the low-energy 
contributions to ah ad > LO in such a way that 91% of 
it comes from the energy region below 1.8 GeV and 
73% is due to the nn channel. Further details can be 
found in Q- 

If isospin were an exact symmetry, we would be 
able to use in ([1]) the spectral functions measured in 
t~ — > X~v decays, where X~ is the (I — 1, J 3 = —1) 
isotopic partner of the X° state. We can define an 
isotopic analogue of the ratio R^ (s) as follows (this 
quantity is related to the usual spectral function 
b Y Rf_(s) = 3v x -(s)): 



Rf.(s) 



B x - 1 dNz 



6 Ka 



B e N x ds 



(2) 



1- 



1 



2s 



In Eq. ([2]), (l/N x )dN x /ds is the normalized invari- 
ant mass spectrum of the hadronic final state, and 
B x - denotes the branching fraction of t — » X~(-f)is T 
(throughout this paper, final state photon radiation 
is implied for r branching fractions). For numerical 
purposes 0], we use for the r lepton mass the value 
m T = (1776.84 ± 0.17) MeV 0, and for the CKM 
matrix element \V ud \ = 0.97418 ± 0.00019 11], which 
assumes CKM unitarity. For the electron branc hing 
fraction we use B e = (17.818 ±0.032)%, obtained [l2| 
supposing lepton universality. 

In the presence of IB effects, the spectral func- 
tion @ in t decays must be corrected by the factor 
Rib(s)/ S E wi in order to be used into the dispersion 
integral ([I]). Therefore, it becomes convenient to in- 
troduce the shift in the dispersion integral: 



Aa h , 



[x- 



3^2 



ds^-Rf. (s) 



4 m 2 



Rm(s) 



(3) 



produced by the IB corrections. The short-distance 
electroweak radiative effects encoded in Sew, which 
includes the re-summation of terms of 0(a n ln n (m z )) 
and of 0(aa™ In"[mz]), lead to the correction S^w = 
1.0235 ± 0.0003 0, EB^lll; the quoted uncertainty 
is attributed to neglected corrections of 0(aa s /Tt) 
15| . This term provides the largest of IB effects in 



-.had.LO 



[X~,t], as it can be seen in Table 1. The re- 
maining IB effects included in Ri B {s) are discussed 
below. 



Hereafter we focus on the X 



channel of 



r lepton decays. Beyond its rather large contribution 
to a^ ad,LO , the precision attained in the measurement 
of the muon anomalous magnetic moment requires 
that the tttt contribution be evaluated below the 1% 
accuracy, making crucial the reliable computation of 
IB corrections [3J, |j] . The s-dependent IB correction 
introduced in (|3]) is defined as: 



Rm(s) 



FSR(s) ffi(s) 
G EM (s) Pi(s) 



F (s) 



F_(s) 



(4) 



The subscripts i = 0, — refer to the electric charge 
of the 2-7T system produced in e + e" annihilation and 
in t~ lepton decays, respectively. Each of the fac- 
tors in i?i B (s) becomes unity in the limit of isospin 
symmetry, thus also i?i B (s) = 1 in this limit. 

In the Standard Model of quarks and leptons in- 
teractions, isospin symmetry is broken by the mass 
difference of u and d quarks, and by the effects of elec- 
tromagnetic interactions. At the hadron level, the IB 
effects introduce some model dependence: hadronic 
matrix elements that are related by isospin symme- 
try, get modified in the presence of IB effects by pho- 
tonic interactions and by the mass and width splitting 
of hadrons involved. Therefore, the usual procedure 
to test isospin symmetry predictions consist in com- 
paring 'bare' hadronic matrix elements obtained from 
experimental data by removing the effects of IB cor- 
rections. 

In the following we consider each of the energy- 
dependent factors that enter in Rm{s) and quantify 
their effects in Aa^ ad,LO [tttt,t]. The correspond- 
ing corrections induced in the branching fraction of 
t — > 7T7W, which is an independent test of these IB 
corrections, can be found in 0, 0j- 

2.1 FSR and phase space corrections 

The final state photonic corrections to e + e~ — > 
7r + 7r~, FSR(s), and the ratio of pion velocities 
Po{s)//3-(s), are the best known corrections to be 
considered in Eq. Q. The FSR correction is com- 
puted using scalar QED and its expression is known 
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Table 1. Contributions to a^ ad,LO [7T7t,t] from the IB corrections discussed in section 2 and ref. [7(. The 
twofold corrections in the second column correspond to results obtained using the GS [3] and KS [3] 



parametrization of pion form factors, respectively. 


For comparison, the last column, denoted 


as DEHZ03, 


contains the results of Ref. [4} . 












Aa h u ^ LO \7rn,r] (IO- 


J 




Source 


CS model 




KS model 


DEHZ03 


Sew 




-12.21±0.15 




-12.1 + 0.3 


Gem 




-1.92±0.90 




-1.0 


FSR 




+4.67±0.47 




+4.5 


p-uu interference 


+2.80±0.19 




+2.80±0.15 


+3.5 + 0.6 


m 7I ± —m^o effect on a 




-7.88 




-7.0 


m 7T ± —m n o effect on F p 


+4.09 




+4.02 


+4.2 


m p ± - m p Q aro 


o.2o±° o: ;; 




11+ 019 

U - LX -0.11 


0.0 + 2.0 


7T7T7, electrom. decays 


-5.91±0.59 




-6.39 + 0.64 


-1.4+1.2 




-16.07±1.22 




-16.70 + 1.23 




Total 




-16.07±1.85 




-9.3 + 2.4 



analytically [16j. From Figure 1 we observe that 
the effects of these two corrections are important close 
to threshold and they vanish rapidly for increasing 
values of s. The phase-space factor is very accurate 
as it depends only on the pion masses. Instead, we 
have attributed a ±10% uncertainty (see Table 1) to 
the contribution of FSR in Aa^ ad,LO [7T7t,t] to account 
for possible deviations from scalar QED. As it has 
been pointed out in @], KLOE [I3| and BABAR [lS | 
measurements of 7r + 7r _ 7(7) in electron-positron col- 
lisions support the validity of this hypothesis within 
the uncertainties quoted above. 

2.2 Long-distance correction 



The definition of the long-distance photonic cor- 
rection Gem to the photon-inclusive hadronic spec- 
trum in r — > ttttv decay can be found elsewhere 
19, The virtual + real soft-photon corrections 
(which gives an infrared-convergent result and that 
we have named mo del- independent corrections in pre- 
vious works [1, of Refs. 0, EiJ are very simi- 
lar numerically, despite the different pion-form factors 
used in both cases (resonance chiral model 2l| and 
vector meson dominance model, respectively). 
This is an expected behavior since G E m is defined 
from the ratio of radiatively-corrected and tree-level 
2ir spectra 0, , thus the form factor dependences 
largely cancel. 

The main difference between the calculations of 
Refs. 



tributions in previous works) of real photon emission, 
and it can be traced back to the model-dependent 
contribution to r — > ixixv^j involving the punr vertex 

20]. In practice, most of the experiments remove 
from their 7r~7r° spectrum, the events associated to 
the decay chain r~ — > 7r~cj(— > 7r°7)^, leaving the in- 
terferences of this with other r lepton radiative am- 
plitudes in the their 7T7t invariant mass distributions 

7]. In order to remain consistent, we also removed 
from our G E m correction the square of the radiative 
amplitude 
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19j and [20( stems from the regular part (which 



is infrared finite, and we call model-dependent con- 



10 1 

s (GeV 2 ) 

Fig. 1. Energy-dependent IB corrections con- 
tained in _Rib(s), Eq. @. 

involving the pun vertex Q. The resulting long- 
distance correction G EM gets closer to the one re- 
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ported in ref. [191 ] ■ The inverse of G EM is plotted 
as a solid line in Figure 1 and its effect in the disper- 
sion integral ((3]) is shown in the second row of Table 
1. We have taken the difference between the effects 
of our Gem in Aa^ ad,LO [7T7t,t] and the one of refer- 
ence [l9| (third column in Table 1) as an estimate 
of the uncertainty associated to model-dependence of 
the long-distance correction. 

2.3 IB effects in pion form factors 

The last IB correction factor in Eq. (j4]), the ratio 
of the electromagnetic to the weak pion form factors, 
involves two sources of IB: (a) a term that mixes the 
1 = 1 and components of the electromagnetic cur- 
rent which is driven by the p — ui mixing and, (b) the 
mass and width difference of p vector mesons which 
affect only the 1 = 1 component of the form factors. 
We discuss this contribution in more detail since it 
represents the main change in Aa^ ad,LO [7r7r,r] with 
respect to previous evaluations of IB corrections. 

Under the above considerations, the pion form fac- 
tors can be written as 0, Q 



F (s) 
F_(s) 



/„-(*). 



l + 6 



ml,-s-im u T u (s) 



,(5) 
(6) 



where the complex parameter 5 pu represents the 
strength of the p — uj mixing, and f p a,-{s) denote 1 " 
the 1 = 1 parts of the pion form factors which are 
dominated, below y/s < 1 GeV, by the p(770) vector 
meson. 

There arc different paramctrizations of the form 
factors fpO,-(s) in the literature which are inspired 
by different models 23j of the p meson propagator. 
However, one would expect that their ratio in Eq. ((4]) 
is relatively less sensitive to a particular model. Just 
for comparison, we adopt two commonly used phe- 
nomenological formulae: the Gounaris-Sakurai (GS) 
24 1 and the Kuhn-Santamaria (KS) [i^ parametriza- 
tions. Consequently, the corrections induced in 
Aa^ ad,LO [7T7r,r] by the IB parameters in the pion 
form factors are quoted as two separate values in the 
second column of Table 1 . 

In the following we discuss the different sources of 
IB in formulae (5) and (6): 

• Strength of p — us mixing: to obtain 6 pu , 
we have fitted Q Eq. (5) with the GS and KS 
parametrizations to the e + e~ data in the full energy 
range available and we have included the effects of 
higher 1=1 resonances in F (s). This approach dif- 
fers from other recent determinations of the p—co mix- 



ing strength which obtain 5 pu from fits to e + e~ data 
below 1 GeV for a wider class of F (s) models [ill ]. 
As a result of our fits reported in [7J, we get for the 
strength and phases of the p — ui mixing parameter: 
\8^\ = (2.00±0.06)xlCT 3 , arg(<^ s ) = (11.6±1.8)°, and 
|^ s | = (1.87±0.06)xl0- 3 , arg(^ s ) = (13.2±1.7)°. In 
both cases, GS and KS, we use an energy- dependent 
absorptive part of the p meson propagator given by 
— iy/sT p o,- (s). Contrary to claims raised in a recent 
paper (23J , we do not find a strong model-dependence 
of the 5 puj mixing parameter. 

• Width difference of p ± — p° mesons 
The energy-dependent decay widths of neutral and 
charged p mesons cannot be measured in an inde- 
pendent way with the accuracy required to estimate 
their effects in Eq. (3). Thus, the width difference 
Ar p = T p ± — T p o must be computed from the total 
widths which are defined as a sum over their exclu- 



sive decay channels 26j. A simple counting of decay 
channels of charged and neutral p mesons give [26| 



Ar p = T[p ± ^n ± n ( 1 )}-T[p a ^7: + n-( 1 )} 

-0.08 MeV, (7) 



where the first two terms include the photon inclu- 
sive rates into two pions (7T7r(7 soft ) + irir~f). The last 
numerical term in Eq. (7) accounts for the rather 
small difference of the remaining decay widths [Io| 
(?T7, 777, 

The 2tt photon inclusive decay rates, first line in 
(7), were calculated including the virtual plus real 
photon radiative corrections in [2ri| . We include its 
energy- dependence in the following form: 



Ar 



7i " 7r ('y) " 



T pTTlT 1 



48tt 



[p_ is){1 + S_)-ffi(s)(l + 6 )] 



(8) 

where g p7r7r denotes the pirn coupling and <5_ >0 con- 
tains the effects of photonic radiative corrections with 
real photons of all energies. Eq. ((U gives Ar^^( 7) = 
(—0.76 ±0.08) MeV at yfs = m p , which can be com- 
pared with a previous estimate, Ar[p — > utt(^)] = 
(+0.49 ±0.58) MeV 0], which was obtained by in- 
cluding only the effects of hard real photon emission 
27J . The ±10% uncertainty added to our result for 
the width difference is estimated from the difference 
observed between our predicted branching fraction for 



"7 



261 ] and its measured value [10 1 



1 These form factors are normalized to unity when s = 0. 
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Fig. 2. Comparison of the ratio of / = 1 com- 
ponents of pion form factors: (upper plot) our 
results (ia ]. and (lower plot) the results of ref. 

The dashed, solid and dashed-dotted lines 
corresponds to m p ± —m„o = (+1,0,-1) MeV. 



As it can be seen from a comparison of the second and 
third columns in Table 1, the width difference (which 
we call '7T7T7, electrom. decays') induces the biggest 
change in a h ad - LO compared to results of previous esti- 
mates Q . Just to emphasize the origin of this impor- 
tant change, in Figure 2 we compare the ratio of I = 1 
components of our form factors 0, @| and the ones 
used in previous calculations 0, for three different 
values of the mass difference: m p +—m p o = (+1,0,-1) 
MeV. The clear difference near the p resonance peak, 
produces the large change in the IB effects due to 7r7r7 
electromagnetic decays. 

• Mass difference of p ± —p° mesons. In previ- 
ous analysis of the IB effects on Aa^ ad,LO Q] it was as- 
sumed that the charged and neutral p mesons were de- 



generated, namely 5m p 



'■ m„±— m„o 



: (0+1) MeV. An 



IB effect in the p meson mass difference arises from 
the self-energy contribution generated by the p° — 7 
mixing term, which affects only the neutral p meson 
mass by m p a —m p o 3r(p° — > e + e~)/{2a) = 1.45 
MeV 71. When we remove this IB effect from the 
value m„± — m„a = (—0.4 + 0.9) MeV measured by 
KLOE [28[ from the Dalitz plot analysis of — ► 
Tr+TT-TT , we find 6m p = (1.0 + 0.9) MeV Q. We use 
this mass splitting in our evaluations shown in Ta- 
ble 1. 



W: 1.075 



1 .05 - EM decays 



£= 

•2 1.025 
o 



O 

O 

CO 



1 = 



0.975 - 



t — | — 1 — 1 — 1 — | — 1 — i — 1 — | — 1 — 1 — 1 — | — 1 — 1 — i — | — i — 1 — r 

— Sm(7T-ji ) 

--■ 5m(p~-p° )=+1.0 MeV 

bare 



J I I I I I I I I I I I I I I L 



_L 



_L 



0.95 

0.2 0.4 0.6 0.8 1 1.2 

s (GeV 2 ) 

Fig. 3. Contributions of mass and width differ- 
ences to the ratio of / = 1 components of form 
factors. 

The IB effect produced in a ^ ad,LO [7T7t,t] by the p 
mass difference is shown in Table 1 for the KS and 
GS parametrizations. As it can be observed, this ef- 
fect gives a rather small contribution. 

In Figure 3 we plot separately the IB corrections 
in the ratio of the 1 = 1 components of the pion form 
factors discussed in subsection 2.3, while their com- 
bined effects including the p~ uj mixing term is repre- 
sented with a dashed-dotted line in Figure 1. The IB 
effects that stems from the ratio of pion form factors 
are important close to the resonance peak, making 
a had,LO [t^-,-] particularly sensitive to these correc- 
tions. The effect of the pion mass difference that 
affects the p meson decay widths is shown in the 6th. 
row of Table 1 and with a solid line in Figure 3. In- 
terestingly, the effects of the pion mass difference and 
of the 77777 electromagnetic decays partly cancel each 
other as it can be observed from Table 1 and Figure 3. 

3 Conclusions 

Isospin breaking (IB) corrections are of great rel- 
evance to improve the accuracy and gain confidence 
on the Standard Model prediction of the leading order 
hadronic contribution a h ad > LO ^ ^ ne mu0 n anomalous 
magnetic moment. These corrections are also impor- 
tant in view of future and more precise measurements 
of the muon anomalous magnetic moment We 
have presented in this work, a summary of some re- 
cent results about these IB corrections. 

As it can be concluded from our results summa- 
rized in Table 1, the new IB corrections produce the 
change AaJ; ad < LO [7T7t,t] = (-16.07±1.85)xl0- 10 which 
is larger by —6.8 x 10~ 10 units when compared to re- 
sults used previously in Q. The main change in the 
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new corrections is due to the effect of the p ± —p° width 
difference 26[ , which quantifies an important IB cor- 
rection near the resonance of the irir system. We have 
calculated, in two commonly used phenomenological 
models 24|, |25[ , the effects of IB corrections that are 
important around the p resonance region, and have 
found that the model-dependence of pion form factors 
is not very important. 

The new IB corrections get closer the results of 
Qhad,LO Dasec { on e + e~ and r lepton data (see 
0j ) ■ These corrections also affect the prediction of 
the r — > TTTrv branching fraction obtained from e + e~ 
data via the isospin symmetry. As it was discussed in 
@,0|; the large discrepancy observed in previous com- 
parisons of CVC predictions and direct measurements 



of this observable Q is also reduced to an acceptable 
level after the new IB corrections are applied. It is 
very appealing that the new IB corrections reduce 
simultaneously the different manifestations of the so- 
called e + e" vs. t lepton discrepancy. 

/ would like to thank Changzheng and his col- 
leagues for their invitation to this workshop and 
very kind hospitality. The results presented here are 
a summary of the work done with M. Davier, A. 
Hocker, B. Malaescu, X. Mo, G. Toledo Sanchez, P. 
Wang, C.Z. Yuan and Z. Zhang. I am very grate- 
ful to them for an enjoyable collaboration and very 
clarifying discussions. 
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